Creating training sequences for space-time diversity arrangements

ABSTRACT

A training sequence is created for space-time diversity arrangement, having any training sequence length, while limiting the training sequence to a standard constellation. Given a number of channel unknowns that need to be estimated, L, a training sequence can be creates that yields minimum means squared estimation error for lengths N t =kN PRUS +L−1, for any positive integer k≧1, where N PRUS  is a selected perfect roots-of-unity sequence (PRUS) of length N. The training sequence is created by concatenating k of the length N perfect roots-of-unity sequences, followed by L−1 initial symbols of that same PRUS. Good training sequences can be created for lengths N t  that cannot be obtained through the above method by concatenating a requisite number of symbols found through an exhaustive search.

RELATED APPLICATIONS

This Application is a Continuation of U.S. patent application Ser. No.11/488,342 filed Jul. 18, 2006. This Application is also related to U.S.patent application Ser. No. 09/957,293 filed Sep. 20, 2001. Thisinvention claims priority from provisional application No. 60/282,647,filed Apr. 9, 2001.

BACKGROUND OF THE INVENTION

This relates to space-time coding, and more particularly, to channelestimation in space-time coding arrangements.

Space-Time coding (STC) is a powerful wireless transmission technologythat enables joint optimized designs of modulation, coding, and transmitdiversity modules on wireless links. A key feature of STC is thatchannel knowledge is not required at the transmitter. While severalnon-coherent STC schemes have been invented that also do not requirechannel information at the receiver, they suffer performance penaltiesrelative to coherent techniques. Such non-coherent techniques aretherefore more suitable for rapidly fading channels that experiencesignificant variation with the transmission block. However, for quasistatic or slowly varying fading channels, training-based channelestimation at the receiver is commonly employed, because it offersbetter performance.

For single transmit antenna situations, it is known that a trainingsequence can be constructed that achieves a channel estimation withminimum mean squared error (optimal sequences) by selecting symbols froman N^(th) root-of-unit alphabet of symbols

${\mathbb{e}}^{\frac{{\mathbb{i}2\pi}\; k}{N}},{k = 0},1,2,\ldots\mspace{14mu},\left( {N - 1} \right),$when the alphabet size N is not constrained. Such sequences are thePerfect Roots-of-Unity Sequences (PRUS) that have been proposed in theliterature, for example, by W. H. Mow, “Sequence Design for SpreadSpectrum,” The Chinese University Press, Chinese University of HongKong, 1995. The training sequence length, N_(t), determines the smallestpossible alphabet size. Indeed, it has been shown that for any givenlength N_(t), there exists a PRUS with alphabet size N=2 N_(t), and thatfor some values of N_(t) smaller alphabet sizes are possible. It followsthat a PRUS of a predetermined length might employ a constellation thatis other than a “standard” constellation, where a “standard”constellation is one that has a power of 2 number of symbols. Binaryphase shift keying (BPSK), quadrature phase shift keying (QPSK), and8-point phase shift keying (8-PSK) are examples of a standardconstellation. Most, if not all, STC systems employ standardconstellations for the transmission of information.

Another known approach for creating training sequences constrains thetraining sequence symbols to a specific (standard) constellation,typically, BPSK, QPSK, or 38-PSK in order that the transmitter andreceiver implementations would be simpler (a single mapper in thetransmitter and an inverse mapper in the receiver—rather than two). Insuch a case, however, optimal sequences do not exist for all traininglengths N_(t). Instead, exhaustive searches must be carried out toidentify sub-optimal sequences according to some performance criteria.Alas, such searches may be computationally prohibitive. For example, inthe third generation TDMA proposal that is considered by the industry,8-PSK constellation symbols are transmitted in a block that includes 116information symbols, and 26 training symbols (N_(t)=26). No optimaltraining sequence exists for this value of N_(t) and constellation sizeand number of channel taps to estimate, L.

When, for example, two transmit antennas are employed, a trainingsequence is needed for each antenna, and ideally, the sequences shouldbe uncorrelated. One known way to arrive at such sequences is through anexhaustive search in the sequences space. This space can be quite large.For example, when employing two antennas, and a training sequence of 26symbols for each antenna, this space contains 8^(2×26) sequences. Forcurrent computational technology, this is a prohibitively large spacefor exhaustive searching. Reducing the constellation of the trainingsequence to BPSK (from 8-PSK) reduces the search to 2^(2×26) sequences,but that is still quite prohibitively large; and the reduction to a BPSKsequence would increase the lowest achievable mean squared error.Moreover, once the two uncorrelated sequences are found, a generator isnecessary for each of the sequences, resulting in an arrangement (for atwo antenna case) as shown in FIG. 1, which includes transmitter 10 thatincludes information encoder 13 that feeds constellation mapper 14 thatdrives antennas 11 and 12 via switches 15 and 16. To provide trainingsequences, transmitter 10 includes sequence generator 5 followed byconstellation mapper 6 that feeds antenna 11 via switch 15, and sequencegenerator 7.

SUMMARY OF THE INVENTION

An advance in the art is achieved with an approach that recognizes thata training sequence can be created for any training sequence length,while limiting the training sequence to a standard constellation.Specifically, given a number of channel unknowns that need to beestimated, L, a training sequence can be creates that yields minimummeans squared estimation error for lengths N_(t)=kN_(PRUS)+L−1, for anypositive integer k≧1, where N_(PRUS) is a selected perfectroots-of-unity sequence (PRUS) of length N, that is smaller than N_(t).The training sequence is created by concatenating k of the length Nperfect roots-of-unity sequences, followed by L−1 initial symbols ofthat same PRUS. Alternatively, k+1 instances of the selected PRUS can beconcatenated, followed by a dropping of symbols in excess of N_(t). Goodtraining sequences can be created for lengths N_(t) that cannot beobtained through the above method by concatenating a requisite number ofsymbols found through an exhaustive search. That is, such lengths can beexpressed by N_(t)=kN_(PRUS)+L−1+M , where k and N_(PRUS) are chosen toyield a minimum M.

BRIEF DESCRIPTION OF THE DRAWING

FIG. 1 shows a prior art arrangement of a two-antenna transmitter and aone antenna receiver, where training sequences are independentlygenerated for the two transmitting antennas;

FIG. 2 shows an arrangement where the training sequences employ the sameconstellation mapper that is employed in mapping space-time encodedinformation symbols;

FIG. 3 presents a block diagram of an arrangement where a single encodergenerates the training sequences for the two transmitter antennas;

FIG. 4 shows one encoding realization for encoder 9 of FIG. 3;

FIG. 5 shows another encoding realization for encoder 9 of FIG. 3;

FIG. 6 presents a block diagram of an arrangement where a single encodergenerates both the training sequences and the information symbols forthe two transmitter antennas;

FIG. 7 shows one encoding realization for encoder 9 of FIG. 6;

FIG. 8 shows another encoding realization for encoder 9 of FIG. 6;

FIG. 9 shows yet another encoding realization for encoder 9 of FIG. 6;

FIG. 10 shows the constellation of an 8-PSK encoder realization forencoder 9; and

FIG. 11 presents a flow diagram for designing L-perfect sequences.

DETAILED DESCRIPTION

The following mathematical development focuses on a system having twotransmit antennas and one receive antenna. It should be understood,however, that a skilled artisan could easily extend this mathematicaldevelopment to more than two transmit antennas, and to more than onereceive antenna.

FIG. 2 shows an arrangement a transmitter with two transmit antennas 11and 12 that transmits signals s₁ and s₂, respectively, and a receiverwith receive antenna 21, and channels h₁ (from antenna 11 to antenna 21)and h₂ (from antenna 12 to antenna 21) therebetween. Channels h₁ and h₂can be expressed as a finite impulse response (FIR) filter with L taps.Thus, the signal received at antenna 21 at time k, y(k), can beexpressed as

$\begin{matrix}{{y(k)} = {{\sum\limits_{i = 0}^{L - 1}{{h_{1}(i)}{s_{1}\left( {k - i} \right)}}} + {\sum\limits_{i = 0}^{L - 1}{{h_{2}(i)}{s_{2}\left( {k - i} \right)}}} + {z(k)}}} & (1)\end{matrix}$where z(k) is noise, which is assumed to be AWGN (additive whiteGaussian noise).

The inputs sequences s_(i) and s₂ belong to a finite signalsconstellation and can be assumed, without loss of generality, that theyare transmitted in data blocks that consist of N, information symbolsand N_(t) training symbols. If N_(t) training symbols are employed toestimate the L taps of a channel in the case of a single antenna, thenfor a two antenna case such as shown in FIG. 1, one needs to employ2N_(t) training symbols to estimate the 2L unknown coefficients (of h₁and h₂).

When a training sequence of length N_(t) is transmitted, the first Lreceived signals are corrupted by the preceding symbols. Therefore, theuseful portion of the transmitted N_(t) sequence is from L to N_(t).Expressing equation (2) in matrix notation over the useful portion of atransmitted training sequence yields

$\begin{matrix}{{y = {{{Sh} + z} = {{\left\lbrack {{S_{1}\left( {L,N_{t}} \right)}\mspace{14mu}{S_{2}\left( {L,N_{t}} \right)}} \right\rbrack\begin{bmatrix}{h_{1}(L)} \\{h_{2}(L)}\end{bmatrix}} + z}}},} & (2)\end{matrix}$where y and z are vectors with (N_(t)−L+1) elements, S₁(L, N_(t)) andS₂(L, N_(t)) are convolution matrices of dimension (N_(t)−L +1)×L , andh₁(L) and h₂(L) are of dimension L×1; that is,

$\begin{matrix}{{{S_{i}\left( {L,N_{t}} \right)} = \begin{bmatrix}{s_{i}\left( {L - 1} \right)} & {s_{i}\left( {L - 2} \right)} & \ldots & {s_{i}(0)} \\{s_{i}(L)} & {s_{i}\left( {L - 1} \right)} & \ldots & {s_{i}(1)} \\\vdots & \vdots & \ddots & \vdots \\{s_{i}\left( {N_{t} - 1} \right)} & {s_{i}\left( {N_{t} - 2} \right)} & \ldots & {s_{i}\left( {N_{t} - L} \right)}\end{bmatrix}},{and}} & (3) \\{{{h_{i}(L)} = {{\begin{bmatrix}{h_{i}(0)} \\{h_{i}(1)} \\\vdots \\{h_{i}\left( {L - 1} \right)}\end{bmatrix}\mspace{14mu}{for}\mspace{14mu} i} = 1}},2.} & (4)\end{matrix}$

If the convolution matrix is to have at least L rows, N_(t) must be atleast 2L−1. In the context of this disclosure, the S matrix is termedthe “training matrix” and, as indicated above, it is a convolutionmatrix that relates to signals received from solely in response totraining sequence symbols; i.e., not corrupted by signals sent prior tothe sending of the training sequence.

The linear least squared channel estimates, ĥ, assuming S has fullcolumn rank, is

$\begin{matrix}{{\hat{h} = {\begin{bmatrix}{\hat{h}}_{1} \\{\hat{h}}_{2}\end{bmatrix} = {\left( {S^{H}S} \right)^{- 1}S^{H}y}}},} & (5)\end{matrix}$where the (□)^(H) and (□)⁻¹ designate the complex-conjugate transpose(Hermitian) and the inverse, respectively. For zero mean noise, thechannel estimation mean squared error is defined byMSE=E[(h−ĥ)^(H)(h−ĥ)]=2σ² tr((S ^(H) S)⁻¹),  (6)where tr(.) denotes a trace of a matrix. The minimum MSE, MMSE, is equalto

$\begin{matrix}{{{MMSE} = \frac{2\sigma^{2}L}{\left( {N_{t} - L + 1} \right)}},} & (7)\end{matrix}$which is achieved if and only if

$\begin{matrix}{{S^{H}S} = {\begin{bmatrix}{S_{1}^{H}S_{1}} & {S_{2}^{H}S_{1}} \\{S_{1}^{H}S_{2}} & {S_{2}^{H}S_{2}}\end{bmatrix} = {\left( {N_{t} - L + 1} \right)I_{2L}}}} & (8)\end{matrix}$where I_(2L) is the 2L×2L identity matrix. The sequences s₁ and s₂ thatsatisfy equation (8) are optimal sequences. Equation (8) effectivelystates that the optimal sequences have an impulse-like autocorrelationfunction (e.g. S₁ ^(H)S corresponds to the identity matrix, I,multiplied by a scalar) and zero cross-correlations.

A straightforward approach for designing two training sequences oflength N_(t) each is to estimates two L-taps channels (i.e., twochannels having L unknowns each, or a total of 2L unknowns) is to designa single training sequence s of length N′_(t) (N′_(t)=N_(t)+L) thatestimates a single channel with L′=2L taps (i.e., a single channelhaving 2L unknowns). Generalizing, N′_(t)=N_(t)+(n−2)L , where n is thenumber of antennas. One can thus view the received signal asy=S(L′,N _(t)′)h(L′)+z  (9)where S is a convolution matrix of dimension (N_(t)′−L′+1)×L′. Again,for optimality, the imposed requirement is thatS ^(H)(L′,N _(t)′)S(L′,N _(t)′)=(N′ _(t) −L′+1)I _(2L),  (10)and once the sequence s is found, the task is to create the subsequencess₁ and s₂ from the found sequence s. Preferably, the subsequences s₁ ands₂ can be algorithmically generated from sequence s. Conversely, one mayfind subsequences s₁ and s₂ that satisfy the requirements of equation(8) and be such that sequence s can be algorithmically generated. Thispermits the use of a single training signal generator that, through apredetermined algorithm (i.e., coding) develops the subsequences s₁ ands₂. Both approaches lead to embodiment depicted in FIG. 3, whereinformation signals are applied to encoder 13 that generates two streamsof symbols that are applied to constellation mapper 14 via switches 15and 16. Generator 5 creates a training sequence that is applied toencoder 9, and encoder 9 generates the subsequences s₁ and s₂ that areapplied to constellation mapper 14 via switches 15 and 16.

Actually, once we realized that the complexity of the training sequencedetermination problem can be reduced by focusing on the creation of asingle sequence from which a plurality of sequences that meet therequirements of equation (8) can be generated, it became apparent thatthere is no requirement for s to be longer than s₁ and s₂.

FIG. 4 presents one approach for generating optimal subsequences s₁ ands₂ that meet the requirements of equation (8) and that can be generatedfrom a single sequence. In accordance with FIG. 4, generator 5 developsa sequence s of length N_(t)/2, and encoder 9 develops therefrom thesequences s₁=−s|s and s₂=s|s , where the “|” symbol stands forconcatenation; e.g., sequence s₁ comprises sequence −s concatenatedwith, or followed by, sequence s. Thus, during the training sequence,antenna 11 transmits the sequence −s during the first N_(t)/2 timeperiods, and the sequence s during the last N_(t)/2 time periods.Antenna 12 transmits the sequence s during both the first and lastN_(t)/2 time periods.

In response to the training sequences transmitted by antennas 11 and 12,receiving antenna 21 develops the signal vector y (where the elements ofthe vector y are the signals received from antennas 11 and 12).Considering the received signal during the first N_(t)/2 time periods asy₁ and during the last N_(t)/2 time periods as y₂, and employing onlythe useful portion of the signal (that is, the portions not corrupted bysignals that are not part of the training sequence) one gets

$\begin{matrix}{\begin{bmatrix}y_{1} \\y_{2}\end{bmatrix} = {{\begin{bmatrix}{- S} & S \\S & S\end{bmatrix}\begin{bmatrix}h_{1} \\h_{2}\end{bmatrix}} + \begin{bmatrix}z_{1} \\z_{2}\end{bmatrix}}} & (11)\end{matrix}$where S is a convolution matrix of dimension (N_(t)−L+1)×L . Inaccordance with the principles disclosed herein, the FIG. 3 receivermultiplies the received signal in processor 25 with the transposeconjugate matrix S^(H), yielding

$\begin{matrix}{\begin{bmatrix}r_{1} \\r_{2}\end{bmatrix} = {{\begin{bmatrix}{- S^{H}} & S^{H} \\S^{H} & S^{H}\end{bmatrix}\begin{bmatrix}y_{1} \\y_{2}\end{bmatrix}} = {{\begin{bmatrix}{2S^{H}S} & 0 \\0 & {2S^{H}S}\end{bmatrix}\begin{bmatrix}h_{1} \\h_{2}\end{bmatrix}} + \begin{bmatrix}{\overset{\_}{z}}_{1} \\{\overset{\_}{z}}_{2}\end{bmatrix}}}} & (12)\end{matrix}$where

$\begin{matrix}{\begin{bmatrix}{\overset{\_}{z}}_{1} \\{\overset{\_}{z}}_{2}\end{bmatrix} = {{\begin{bmatrix}{- S^{H}} & S^{H} \\S^{H} & S^{H}\end{bmatrix}\begin{bmatrix}z_{1} \\z_{2}\end{bmatrix}}.}} & (13)\end{matrix}$If the sequence s is such that S^(H)S=(N_(t)−L+1)I_(L), then

$\begin{matrix}{\begin{bmatrix}r_{1} \\r_{2}\end{bmatrix} = {{2{\left( {N_{t} - L + 1} \right)\begin{bmatrix}h_{1} \\h_{2}\end{bmatrix}}} + {\left\lbrack \begin{bmatrix}{\overset{\_}{z}}_{1} \\{\overset{\_}{z}}_{2}\end{bmatrix} \right\rbrack.}}} & (14)\end{matrix}$If the noise is white, then the linear processing at the receiver doesnot color it, and the channel transfer functions correspond toh ₁=1/2(N _(t) −L+1)r ₁h ₂=1/2(N _(t) −L+1)r ₂  (15)with a minimum squared error, MSE, that achieves the lower boundexpressed in equation (7); to wit,

$\begin{matrix}{{MSE} = {\frac{2\sigma^{2}L}{\left( {N_{t} - L + 1} \right)}.}} & (16)\end{matrix}$

The above result can be generalized to allow any matrix U to be used toencode the training sequence, s, so that

$\begin{matrix}{\begin{bmatrix}y_{1} \\y_{2}\end{bmatrix} = {{U\begin{bmatrix}h_{1} \\h_{2}\end{bmatrix}} + \begin{bmatrix}z_{1} \\z_{2}\end{bmatrix}}} & (17)\end{matrix}$as long as U^(H)U=2I for a two antennas case, and U^(H)U=KI for an Kantenna case.

Whereas FIG. 4 presents a method for developing sequences s₁ and s₂ oflength N_(t) from a sequence s that is N_(t)/2 symbols long, FIG. 5presents a method for developing sequences s₁ and s₂ of length N_(t)from a sequence s that is 2N_(t) symbols long, which consists of asequence d₁=[s(0) s(1) . . . s(N_(t)−1)] followed by a sequence d₂=[s(0)s(1) . . . s(N_(t)−1)]. In accordance with this approach,s_(t)=d₁|{tilde over (d)}*₂ and s₂=d₂|{tilde over (d)}*₁. The sequence{tilde over (d)}₁ corresponds to the sequence d₁ with its elements inreverse order. The symbol {tilde over (d)}*₁ operation corresponds tothe sequence d₁ with its elements in reverse order and converted totheir respective complex conjugates.

The FIG. 5 encoding is very similar to the encoding scheme disclosed byAlamouti in U.S. Pat. No. 6,185,258, issued Feb. 6, 2001, except that(a) the Alamouti scheme is symbols-centric whereas the FIG. 5 encodingis sequence-centric, and (b) the Alamouti scheme does not have theconcept of a reverse order of a sequence (e.g., {tilde over (d)}*₁). Seealso E. Lindskog and A. Paulraj titled “A Transmit Diversity Scheme forChannels With Intersymbol Interference,” ICC, 1:307-311, 2000. Anencoder 9 that is created for developing training sequences s₁ and s₂ inaccordance with FIG. 5, can be constructed with a control terminal thatis set to 1 during transmission of information and set to another value(e.g., 0, or to N_(t) to indicate the length of the generated block)during transmission of the training sequence, leading to the simplifiedtransmitter realization shown in FIG. 6. More importantly, such anarrangement leads to a simplified receiver because essentially the samedecoder is used for both the information signals and the trainingsignals.

With a signal arrangement as shown in FIG. 5, the signal captured atantenna 21 of receiver 20 is

$\begin{matrix}{\begin{bmatrix}y_{1} \\y_{2}\end{bmatrix} = {{\begin{bmatrix}{- {\overset{\sim}{D}}_{2}^{*}} & {\overset{\sim}{D}}_{1}^{*} \\D_{1} & D_{2}\end{bmatrix}\begin{bmatrix}h_{1} \\h_{2}\end{bmatrix}} + \left\lbrack \begin{bmatrix}z_{1} \\z_{2}\end{bmatrix} \right\rbrack}} & (18)\end{matrix}$where the matrices D_(i) and {tilde over (D)}_(i) (for i=1,2) areconvolution matrices for d₁ and {tilde over (d)}₁, respectively, ofdimension (N_(t)−L+1)×L . Recalling from equation (8) that MMSE isachieved if and only if D^(H)D has zeros off the diagonal; i.e.,−{tilde over (D)} ₁ ^(T) {tilde over (D)}* ₂+(D* ₂)^(T) D ₁=0  (19)and−{tilde over (D)} ₂ ^(T) {tilde over (D)}* ₁+(D* ₁)^(T) D ₂=0,  (20)and identity matrices on the diagonal; i.e.,{tilde over (D)} ₂ ^(T) {tilde over (D)}* ₂+(D* ₁)^(T) D ₁=2(N _(t)−L+1)I _(L)  (21)and{tilde over (D)} ₁ ^(T) {tilde over (D)}* ₁+(D* ₂)^(T) D ₂=2(N _(t)−L+1)I _(L).  (22)

Various arrangements that interrelate sequences d₁ and d₂ can be foundthat meet the above requirement. By way of example (and not by way oflimitation), a number of simple choices satisfy these conditions follow.

-   (1) (D*₁)^(T) D₁=(N_(t)−L+1)I_(L), {tilde over (D)}₁=D₁, and D₂=D₁.    To show that equation (21) holds, one may note that {tilde over    (D)}₂ ^(T){tilde over (D)}*₂ (the first term in the equation)    becomes D₁ ^(T)D*₁, but if (D*₁)^(T) D₁ is a diagonal matrix then so    is {tilde over (D)}₂ ^(T){tilde over (D)}*₂. Thus, according to this    training sequence embodiment, one needs to only identify a sequence    d₁ that is symmetric about its center, with an impulse-like    autocorrelation function, and set d₂ equal to d₁. This is shown in    FIG. 7.-   (2) (D*₁)^(T) D₁=(N_(t)−L+1)I_(L), and {tilde over (D)}₂=D₁. To show    that equation (21) holds, one may note that the {tilde over (D)}₂    ^(T){tilde over (D)}*₂ first term in the equation also becomes D₁    ^(T)D*₁. Thus, according to this training sequence embodiment, one    needs to only identify a sequence d₁ with an impulse-like    autocorrelation function, and set d₂ equal to {tilde over (d)}₁ .    This is shown in FIG. 8.-   (3) (D*₁)^(T) D₁=(N_(t)−L+1)I_(L), and {tilde over (D)}*₂=D₁. To    show that equation (21) holds, one may note that the {tilde over    (D)}₂ ^(T){tilde over (D)}*₂ first term in the equation becomes    (D*₁)^(T) D₁. Thus, according to this training sequence embodiment,    one needs to only identify a sequence d₁ with an impulse-like    autocorrelation function, and set d₂ equal to {tilde over (d)}*₁.    This is shown in FIG. 9.    Training Sequences Employing Trellis Coding

Consider a trellis code with m memory elements and outputs from aconstellation of size C, over a single channel with memory2^(m)C^((L−1))−1. To perform joint equalization and decoding one needs aproduct trellis with 2^(m)C^((L−1)) states. For a space-time trelliscode with m memory elements, n transmit antennas and one receiveantenna, over a channel with memory (L−1), one needs a product trelliswith 2^(m)C^(n(L−1)) .

The receiver can incorporate the space-time trellis code structure inthe channel model to create an equivalent single-input, single outputchannel, h_(eq), of length m+L. The trellis, in such a case, involvesC^((m+L−1)) states. The approach disclosed herein uses a single trainingsequence at the input of the space-time trellis encoder to directlyestimate h_(eq) used by the joint space-time equalizer/decoder. Thechannel h_(eq) that incorporates the space-time code structure typicallyhas a longer memory than the channel h₁ and h₂ (in a system where thereare two transmitting antennas and one receiving antenna).

To illustrate, assume an encoder 30 as depicted in FIG. 10 that employsan 8-PSK constellation of symbols to encode data from a trainingsequence generator into a sequence s of symbols taken from the sete^(i2πp) ^(k) ^(/8), p_(k)=0,1,2, . . . ,7 , where the trainingsequences s₁ and s₂ are algorithmically derived within encoder 30 fromsequence s. Specifically, assume that s₁(k)=s(k), and thats₂(k)=(−1)^(p) ^(k−1) s(k−1), which means that s₂(k)=s(k−1) when s(k−1)is an even member of the constellation (e^(i0), e^(iπ/2), e^(iπ), ande^(i3π/2)), and s₂(k)=−s(k−1) when s(k−1) is an odd member of theconstellation.

With such an arrangement, the received signal at time k can be expressedas

$\begin{matrix}\begin{matrix}{{y(k)} = {{\sum\limits_{i = 0}^{L - 1}\;{{h_{1}(i)}{s\left( {k - i} \right)}}} + {\sum\limits_{i = 0}^{L - 1}\;{{h_{2}\left( {i,k} \right)}\left( {- 1} \right)^{p_{k - i - 1}}{s\left( {k - i - 1} \right)}}} + {z(k)}}} \\{{= {{\sum\limits_{i = 0}^{L}\;{{h_{eq}(i)}{s\left( {k - i} \right)}}} + {z(k)}}},}\end{matrix} & (23)\end{matrix}$where

$\begin{matrix}{{h_{eq}\left( {i,k} \right)} = \left\{ \begin{matrix}{h_{1}(0)} & {{{for}\mspace{14mu} i} = 0} \\{{h_{1}(i)} + {\left( {- 1} \right)^{p_{k - i}}{h_{2}\left( {i - 1} \right)}}} & {{{for}\mspace{14mu} 0} < i < L} \\{\left( {- 1} \right)^{p_{k - L}}{h_{2}\left( {L - 1} \right)}} & {{{for}\mspace{14mu} i} = {L.}}\end{matrix} \right.} & (24)\end{matrix}$A block of received signals (corresponding to the useful portion of thetraining sequence block) can be expressed in matrix form byy=Sh _(eq) +z  (25)where

$\begin{matrix}{S = {\begin{bmatrix}{s(L)} & {s\left( {L - 1} \right)} & \ldots & {s(0)} \\{s\left( {L + 1} \right)} & {s(L)} & \ldots & {s(1)} \\\vdots & \vdots & \; & \vdots \\{s\left( {N_{t} - 1} \right)} & {s\left( {N_{t} - 2} \right)} & \ldots & {s\left( {N_{t} - L - 1} \right)}\end{bmatrix}\begin{bmatrix}{h_{eq}\left( {0,L} \right)} \\{h_{eq}\left( {1,{L + 1}} \right)} \\\vdots \\{h_{eq}\left( {L,{N_{t} - 1}} \right)}\end{bmatrix}}} & (26)\end{matrix}$and following the principles disclosed above, it can be realized thatwhen the training sequence is properly selected so that S^(H)S is adiagonal matrix, i.e., S^(H)S=(N_(t)−L)I_(L+1), an estimate of h_(eq),that is, ĥ_(eq), is obtained from

$\begin{matrix}{{\hat{h}}_{eq} = {\frac{S^{H}y}{\left( {N_{t} - L} \right)}.}} & (27)\end{matrix}$If the training sequence were to comprise only the even constellationsymbols, e^(i2πk/8), k=0,2,4,6 , per equation (24), the elements ofĥ_(eq) would correspond toh _(eq) ^(even) =[h ₁(0), h ₁(1)+h ₂(0), h ₁(2)+h ₂(1), . . . h ₁(L−1)+h₂(L−2), h ₂(L−1)].  (28)If the training sequence were to comprise only the odd constellationsymbols, e^(i2πk/8), k=1,3,5,7, the elements of ĥ_(eq) would correspondtoh _(eq) ^(odd) =[h ₁(0), h ₁(1)−h ₂(0), h ₁(2)−h ₂(1), . . . h ₁(L−1)−h₂(L−2),−h ₂(L−1)].  (29)If the training sequence were to comprise a segment of only evenconstellation symbols followed by only odd constellation symbols (orvice versa), then channel estimator 22 within receiver 20 can determinethe h_(eq) ^(even) coefficients from the segment that transmitted onlythe even constellation symbols, and can determine the h_(eq) ^(odd)coefficients from the segment that transmitted only the evenconstellation symbols. Once both h_(eq) ^(even) and h_(eq) ^(odd) areknown, estimator 22 can obtain the coefficients of h₁ from

$\begin{matrix}{\left\lbrack {{h_{1}(0)},{h_{1}(1)},{h_{1}(2)},{\ldots\mspace{14mu}{h_{1}\left( {L - 1} \right)}}} \right\rbrack = \frac{h_{eq}^{\;{even}} + h_{eq}^{\;{odd}}}{2}} & (30)\end{matrix}$and the coefficients of h₂ from

$\begin{matrix}{\left\lbrack {{h_{2}(0)},{h_{2}(1)},{h_{2}(2)},{\ldots\mspace{14mu}{h_{2}\left( {L - 1} \right)}}} \right\rbrack = {\frac{h_{eq}^{\;{even}} - h_{eq}^{\;{odd}}}{2}.}} & (31)\end{matrix}$What remains, then, is to create a single training sequence s of lengthN_(t) where one half of it (the s_(even) portion) consists of only evenconstellation symbols (even sub-constellation), and another half of it(the S_(odd) portion) consists of only odd constellation symbols (oddsub-constellation). The sequences s₁ and s₂ of length N_(t) are derivedfrom the sequence s by means of the 8-PSK space-time trellis encoder.The sequences s₁ and s₂ must also meet the requirements of equation (8).Once s_(even) is found, s_(odd) can simply bes_(odd)=αs_(even), where α=e^(iπk/4) for any k=1,3,5,7.  (32)Therefore, the search for sequence s is reduced from a search in the of8^(N) ^(t) to a search for s_(even) in the space 4^((N) ^(t) ^(/2)) suchthat, when concatenated with s_(odd) that is computed from s_(even) asspecified in equation (32), yields a sequence s that has anautocorrelation function that is, or is close to being, impulse-like.

For a training sequence of length N_(t)=26, with an 8-PSK space-timetrellis encoder, we have identified the 12 training sequences specifiedin Table 1 below.

TABLE 1 sequence # α Se 1 exp(i5π/4) −1 1 1 1 1 −1 −i −1 1 1 −1 1 1 2exp(i3π/4) 1 1 −1 1 i i 1 −i i −1 −1 −1 1 3 exp(iπ/4) 1 −1 −1 −i i −i 11 1 −i −1 1 1 4 exp(iπ/4) 1 −1 −1 −i 1 −1 1 −i −i −i −1 1 1 5 exp(iπ/4)1 i 1 1 i −1 −1 i 1 −1 1 i 1 6 exp(i3π/4) 1 i 1 i −1 −1 1 −1 −1 i 1 i 17 exp(i7π/4) 1 −i 1 1 −i −1 −1 −i 1 −1 1 −i 1 8 exp(i5π/4) 1 −i 1 −1 i 1−1 −i 1 1 1 −i 1 9 exp(i3π/4) −1 1 1 1 −1 −1 −i −1 −1 1 −1 1 1 10exp(i7π/4) −1 i −1 −i 1 −i i i 1 i 1 −i 1 11 exp(iπ/4) −1 −i −1 i 1 i −i−i 1 −i 1 i 1 12 exp(i3π/4) −1 −i −1 i −1 i −i −i −1 −i 1 i 1Construction of Training Sequence

While the above-disclosed materials provide a very significantimprovement over the prior art, there is still the requirement ofselecting a sequence s₁ with an impulse-like autocorrelation function.The following discloses one approach for identifying such a sequencewithout having to perform an exhaustive search.

A root-of-unity sequence with alphabet size N has complex roots of unityelements of the form

${\mathbb{e}}^{\frac{{\mathbb{i}}\; 2\;\pi\; k}{N}},{k = 1},2,{\ldots\mspace{14mu}{\left( {N - 1} \right).}}$As indicated above, the prior art has shown that perfect roots-of-unitysequences (PRUS) can be found for any training sequence of length N_(t),as long as no constraint is imposed on the value of N. As also indicatedabove, however, it is considered disadvantageous to not limit N to apower of 2. Table 2 presents the number of PRUSs that were found toexist (through exhaustive search) for different sequence lengths whenthe N is restricted to 2 (BPSK), 4 (QPSK), or 8 (8-PSK). Cell entries inTable 2 with “-” indicate that sequence does not exist, and blank cellsindicate that an exhaustive search for a sequence was not performed.

TABLE 2 N_(t)=  2 3  4 5 6 7  8 9 10 11 12 13 14 15  16 17 18 BPSK — — 8 — — — — — — — — — — — — — — QPSK  8 —  32 — — — 128 — — — — — — —6144 — — 8-PSK 16 — 128 — — — 512 — — — — —

A sequence s of length N_(t) is called L-perfect if the correspondingtraining matrix S of dimension (N_(t)−L+1)×L satisfies equation (8).Thus, an L-perfect sequence of length N_(t) is optimal for a channelwith L taps. It can be shown that the length N_(t) of an L-perfectsequence from a 2^(p)-alphabet can only be equal to

$\begin{matrix}{{N_{t} = \begin{Bmatrix}{2\left( {L + i} \right)} & {{{for}\mspace{14mu} L} = {odd}} \\{{2\left( {L + i} \right)} - 1} & {{{for}\mspace{14mu} L} = {even}}\end{Bmatrix}},{{{for}\mspace{14mu} i} = 0},1,\ldots\mspace{14mu},} & (33)\end{matrix}$which is a necessary, but not sufficient, condition for L-perfectsequences of length N_(t). Table 3 shows the minimum necessary N_(t) forL=2,3, . . . 10, the size of the corresponding matrix S, and the resultsof an exhaustive search for L-perfect sequences (indicating the numberof such sequences that were found). Cell entries marked “x” indicatethat sequences exist, but number of such sequences it is not known.

TABLE 3 L 2 3 4 5 6 7 8 9 10 N_(t) 3 6 7 10 11 14 15 18 19 S 2 × 2 4 × 34 × 4 6 × 5 6 × 6 8 × 7 8 × 8 10 × 10 × 9 10 BPSK 4 8 8 — — — — — — QPSK16 64 64 — — 128 x 8-PSK 64 512 512 x x

It is known that with a PRUS of a given length, N_(PRUS), one canestimate up to L=N_(PRUS) unknowns. It can be shown that a trainingsequence of length N_(t) is also an L-perfect training sequence ifN _(t) =kN _(PRUS) +L−1 and k≧1.  (34)Accordingly, an L-perfect sequence of length kN_(PRUS)+L−1 can beconstructed by selecting an N_(PRUS) sequence, repeating it k times, andcircularly extending it by L−1 symbols. Restated and amplified somewhat,for a selected PRUS of a given N_(PRUS), i.e.,s _(p)(N _(PRUS))=[s _(p)(0) s _(p)(1) . . . s _(p)(N _(PRUS)−1)],  (35)the L-perfect sequence of length kN_(PRUS)+L−1 is created from aconcatenation of k s_(p)(N_(PRUS)) sequences followed by the first L−1symbols of s_(p)(N_(PRUS)), or from a concatenation of the last L−1symbols of s_(p)(N_(PRUS)) followed by k s_(p)(N_(PRUS)) sequences.

FIG. 11 shows the simplicity of this method. In step 100 an N_(PURS)sequence is selected in response to an applied training sequence lengthN_(t) (that is greater than nL−1, where n is the number of transmitterantennas). Control then passes to step 101 where the training sequenceis constructed, as indicated above, from concatenation of k of theselected N_(PRUS) sequence followed by the first L−1 symbols of s_(P)(N_(PRUS)), or from a concatenation of the last L−1 symbols of theselected N_(PRUS) sequence followed by k concatenated sequences.

To illustrate, assume that the number of channel “taps” that need to beestimated, L, is 5, and that a QPSK alphabet is desired to be used. Fromthe above it is known that N_(PRUS) must be equal to or greater than 5,and from Table 2 it is known that the smallest N_(PRUS) that can befound for QSPK that is larger than 5 is N_(PRUS)=8. Employing equation(34) yields

$\begin{matrix}\begin{matrix}{N_{t} = {{kN}_{PRUS} + L - 1}} \\{= {{k \cdot 8} + 5 - 1}} \\{{= 12},20,28,{{\ldots\mspace{14mu}{for}\mspace{14mu} k} = 1},2,{3\mspace{14mu}\ldots}}\end{matrix} & (36)\end{matrix}$

While an L-perfect training sequence cannot be constructed from PRUSsequences for values of N_(t) other than values derived by operation ofequation (34), it is known that, nevertheless, L-perfect sequences mayexist. The only problem is that it may be prohibitively difficult tofind them. However, in accordance with the approach disclosed below,sub-optimal solutions are possible to create quite easily.

If it is given that the training sequence is N_(t) long, one can expressthis length byN _(t) =kN _(PRUS) +L−1+M, where M>0  (37)In accord with our approach, select a value of N_(PRUS) ≧L thatminimizes M, create a sequence of length kN_(PRUS)+L−1 as disclosedabove, and then extend that sequence by adding M symbols. The M addedsymbols can be found by selecting, through an exhaustive search, thesymbols that lead to the lowest estimation MSE. Alternatively, select avalue of N_(PRUS)>L that minimizes M′ in the equation,N _(t) =kN _(PRUS) +L−1−M′, where M′>0  (38)create a sequence of length kN_(PRUS)+L−1 as disclosed above, and thendrop the last (or first) M′ symbols.

The receiver shown in FIG. 2 includes the channel estimator 22, whichtakes the received signal and multiplies it S^(H) as appropriate; seeequation (12), above.

1. A method executed in a hardware module for creating a trainingsequence of symbols, of chosen length N_(t), for estimating L unknownsof a wireless channel, comprising the steps of: (a) if N and k existsuch that N<N_(t) is a power of 2, k is a positive integer, andkN+(L−1)=N_(t), concatenating k instances a perfect roots-of-unitysequence of alphabet N (N_(PRUS) sequence), and adding an L−1 portion ofthe concatenated k instances of said N_(PRUS) sequence; (b) if N and kdo not exist such that N<N_(t) is a power of 2, k is a positive integer,and kN+(L−1)=N_(t) (i) selecting k that minimizes M in kN+(L−1)+M=N_(t)where M is a positive integer, concatenating k instances a perfectroots-of-unity sequence of alphabet N (N_(PRUS) sequence), adding an L−1portion of the concatenated k instances of said N_(PRUS) sequence; andadding M additional symbols; or (ii) selecting k that minimizes M inkN+(L−1)−M=N_(t) where M is a positive integer, concatenating k-1instances a perfect roots-of-unity sequence of alphabet N (N_(PRUS)sequence), adding an N+L−1−M portion of the concatenated k instances ofsaid N_(PRUS) sequence; and (c) transmitting said training sequence oversaid channel.